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ABSTRACT 


Consider the two-factor crossed components of variance model given 
hy 

Y iJ = V + A i + T J + e ij for i = 1 * -- 1 * J" 1 . J - 

The random variables A^, Tj, are normal and independent with means zero 

and variances o^, a*. Approximate confidence intervals are presented and 

evaluated for the ratios p a = a*/(o| + + o|) and P T * a2/(a| * + o|)» 







1. Introduction . 

Consider the two-factor crossed Components-of-Variance model 
given by 

Y ij = A i + T j + c ij j=1.2,...0 . 

The ,Tj» j are independent unobservable random variables and 
A. ~ -M(0,o a 2 ) ^ ~ N(0.a T 2 ) - N(0,a £ 2 ) 

H is an unknown parameter, and the Y^j are observable random 
variables. 

The analysis of variance table for this model is given in Table l. 

Table 1. 


Source 

d.f. 

S.S. 

M.S. 

E.H.S. 

Total 

I J 

i.inT 





Mean 

1 

I J Y 2 





Due to A 

n, - (I-D 


i 2 


* 

"i * J % 2 * \ 

Due to T 

n 2 = (J-l) 


i 2 


4 

°2 * l0 T 2 * \ 

Error 

n 3 = V? 


-v 7 .. 

,) 2 

4 

e 3 = ° c 2 


y 
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The random variables Y, S 1 , S 2 , S 3 are complete sufficient 
statistics for this problem and n j S i 2 /ej~x 2 (n. j )for i=l,2,3 and 
are mutually independent. 

An important problem in applied statistics is to obtain confidence 
Intervals on various functions of the variance components o^ 2 , Oj 2 , and 
%*• 

For some well-known functions exact-size confidence intervals a*e 
available, but in mar.y situations exact-size confidence intervals have 
not as yet been derived. 

For the following important and useful functions exact-size 
confidence intervals are known (exact-size means exact confidence 
coefficient): o £ 2 , o*/o* t o. 2 /a T 2 , o e 2 /(o £ 2 + <>/). ° c 2 /(o e 2 + ° T 2 ), 
o^ 2 /(o c 2 + o^ 2 ), o T 2 /(u e 2 + o^. 2 ). For a discussion of these see Cre.ybiii (1976). 

Very good approximate confidence intervals are known for o ^ 

2 

and Oj 

For approximate confidence intervals on linear combinations of 
2 2 2 

, Oj , 0 £ see Welch (1956). 

In applied work it is often important to determine the ratios of 
2 2 2 

the Individual variances , Oy , and to the total variation 

(a^ 2 + Oj 2 + 0 £ 2 ), i.e., to determine 
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There Is no method available for setting exact 1-a confidence 
Intervals on , py , or p £ . Large sample approximations have 
been proposed by Osborne and Paterson (1952). 

The purpose of this article is to obtain good approximate 1-a 
upper, lower, and two-sided confidence intervals for and Py 
and evaluate how close the approximate confidence coefficients arc to the 
specified 1-a . 

The quantities p^ and py are useful In many applications. 

Graybill (1976) gives an example of a two-factor crossed components of 
variance model without interaction. An investigator wants to examine 
the effects of drivers and automobiles on gas mileage. One type of 
automobile is used and 20 are selected at random from one month's 
production. Ten people are selected at random from a large city and 
are asked to drive each car for two weeks and record the number of miles 
per gallon of gas used for each car for the two-week period. The total 
variation In the number of miles per gallon of gas used for this type of 
automobile is °j( + °f + where is the variation due to the cars, 
o T is the variation due to the drivers, and <r is the variation due to 

T C 

all other noncontrolled factors. The proportion of the variation in the 
number of miles per gallon due to the cars is given by , and the 
proportion of the variation in the number of miles per gallon duo to the 
drivers is given by py. 
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2• Approximate Confidence Intervals . 

In this section approximate confidence intervals on p A and 
pj will be derived. 

First consider a 1-a upper confidence interal on p A 
denoted by i^ < o A where «, A is the lower confidence point. 

It is simpler to find a 1-a upper confidence Interval for 

.. y »i -°3 

A I(o T Z + o c 2 ) 0 2 + "l e 3 

and convert this to a confidence interval for p a by using the fact 



A function g L (Y) of the observable random vector Y_ must be 
determined such that 

P[9 L (Y) < e A ] 

Is close to a specified 1-a for all values of the unknown parameters. 
Since Y, S^, S 2 2 , S 3 2 is a set of complete sufficient statistics 
for this problem, the lower confidence point should be a function of 
these random variables. So a function h L (Y,S^ 2 ,S 2 2 ,S 3 2 ) must be 
determined such that 

P[h L (T.s, 2 ,s 2 z ,s 3 2 ) i e A ] 

Is close to a speclfiea 1-a . 

If any constant k is added to each observation the 

2 ? 2 

parameter is unchanged so it is required that h^Y+k.Sj ) 




6 


remains unchanged. The constant k = - Y Is chosen and thus 
h L (T > S 1 2 ,S 2 2 ,S 3 2 ) = h L (0,S 1 2 ,S 2 2 ,S 3 2 ) = <\ l (S*,S 2 2 ,S 3 Z ), which is a 
function of Sj 2 ,S 2 2 , and S 3 2 only. 

If each observation Y^j is multiplied by a non-zero constant 
c , the parameter Is unchanged so it is required that 
q L ( s l 2 , s 2 2 , s 3 2 ) remai ns unchanged also. The constant chosen is 
1/S 3 2 so 

2 

which is a function of $j 2 /S 3 2 = F-j and S 2 2 /S 3 2 * F 2 only. 

The problem is reduced to determining a function L(Fj,F 2 ) such 

that 

P[L(F r F 2 ) < e A ] 

will be close to a specified confidence coefficient 1-a . 

To determine the function L(Fj,F 2 ) the following five conditions 
are imposed. 

1) The confidence interval must be exact when e 3 = 0. 

2 

When e 3 = 0 it follows that S 3 = 0 with probability one, 

Fj"^ = 0 with probability one, and = 0j/e 2 • S° an exfl ct 1-a 
lower confidence limit on e A is 


.S 2 2 .s 3 2 ) 




/s , 2 




a:n, ,n,j2 


and L(F^,F 2 ) must be equal to 
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In order to satisfy this condition. 

2) When o ^ -*» it is required that the confidence coefficient 
2 2 2 

be exact. When -*» and oj ,o £ are fixed it follows that 

-*• with 02 ,e 3 ^ xed » and 8 ] dominates e A . So the confidence 
Interval on e A should behave like an exact 1-a confidence interval 
on 6j . So L(Fj,F 2 ) should be equal to Sj 2 F a> ”^ 


3) Tne confidence interval is required to be exact when 
0 -*» and I is fixed. 

When J -*» it follows that n 2 -*», n 3 -*» and hence S 2 2 -*• e 2 
In probability and S 3 2 ♦ e 3 in probability; also, 

(S 2 2 + ni S 3 2 ) - (e 2 + ni e 3 ) 

in probability. 

Begin with an exact (1-a) confidence interval for e 1 given 

by 

and subtract S 3 2 from the left side and the "equivalent" value e. 


PI 






from the right side to obtain 


S 4 F •„ 

1 a:n r - 


3 - °1 


Then divide the left and right sides, respectively, by 
(S 2 2 + n ^ S 3 2 ) and the "equivalent" value (e 2 + n^) to obtain 


°1 r a:n,,« 


"r3 4 ' J 

Hence, when J-*» the lower confidence point should behave like 


l(f,,f 2 ) - — r 


4) When I-+® and 0 is fixed the confidence interval is 
required to be exact. 

For I-x» and J fixed it follows that , n 3 *» and 

2 ? 

S-j -*■ e 1 in probability, and S 3 -* e 3 in probability; also, 

- S 3 2 -*• e 1 - e 3 in probability. 

Begin with a 1-a exact confidence interval for e 2 given by 

n . C 2 C 


Add njS 3 to the right side and the "equivalent" value n^e 3 to the 
left side to obtain 

e, + n,©, < S, Z F . „ + n, S 2 

Z 1 3 — Z o:»,n 2 1 3 

Take the inverse of each side, then multiply the right and left sides 

2 p 

respectively, by S-j - S 3 and its "equivalent” value - c 3 






to obtain 
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Hence, when I->«,L(F 1 .F,,) is required to behave like 
F, - 1 

L < F r F 2 > ’ nr- 1 —rj- 


accepted the confidence interval should include "zero" and L(F-j .Fg) 

should be increasing in F-j for fixed values of F 2 . 

H is accepted iff F, < F , so the confidence limit 

o i — a. n-j ,nj 

should be 

L(F,,F 9 ) * 0 when F, < F , 

12 l — a:n.j ,n^ 

UF-j.Fg) monotonic increasing in F 1 for fixed values of ? 2 
when F, > F 

i a:nj,n.j . 


a„ + F-, 

« F i' f i»-rnx 

0 1 2 

but a more general function will be examined, i.e., 


where a^ and b^ which are functions of nj^.n^.n are determined 
by requiring L(F.j,F 2 ) to satisfy conditions 1), 2), 3), 4), 5) 
Substituting condition 1) gives 


( 2 . 2 ) 
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Substituting condition 2) gives 


1 a:n,,® 


(2.3) 


^1 F a:n, ,® F a:n,,n 9 


(2.4) 


From condition 3) Lim a = -1 , Lim b = n, , Lim a, = F ~ „ 

J-xo 0 J-x» 0 1 J-x» 1 a:n l’ 


Lim b, * 1, Lim a, * 0 , Lim b 9 = 0 

J-x» 1 J-x* c J-x» c 

which are satisfied for those values of a^, b-j from equations 
(2.3) and (2.4). 

From condition 4) Lim a = - 1, Lim = n, , Lim a, = 1 , 
I- 0 I- 0 1 I- 1 

Lim b, = F „ , Lim a 0 = 0 , Lim b 0 = 0 . 

I- 1 a: “’ n 2 I- 2 I- 2 


and 4). 
From 


condition 5) L(F 1> F 2 ) = 0 if F^ <_ F o<n n which implies 


a o * a l F 1 * a 2 F l ‘ 0 if F 1 i F o; 


Now substitute F, = F _ into 
1 a.n r n 3 


a„ + a,F n + a,F 


-1 


o 1 a:n, ,n_ 2 a.'n, ,n. 


! 0 , 


"r"3 4 1 * 3 

let a-j be the value given in equation (2.3) and solve for a 2 to 


a 2 ' ( " a o ' F a:n rB . F a:n r n 3 )F a:n r n 3 


obtain 
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Use the simplest values for a and b , namely, a = -1 


and b = n, , and obtain 


“1, a ‘ n l *' 3 ___], 3 


1 o:n,,® a:n,,n 9 2 22 


L(F,,F 2 ) = 0 If F,£F a;n]> „ 3 

To determine the value of b^ we impose condition C). 

2 2 

Condition 6j, When H q :oj =0 vs. >0 is accepted 

the confidence interval is required to be exact. 

2 2 

To examine condition 7) we get H q :o t - 0 vs. H^:o t >0 <~> 

H o : 0 2 = 6 3 VS * H l :0 2 * 0 3 w H o :0 A = ( 0 r e 3 )/IO 3 VS ‘ 

H 1 :e A < (o r e 3 )/le 3 . 

If H q is true then e A = Jo' A 2 /Io c 2 = ]- - lj 

So if H q is accepted at a level we want L(F^ .F^) to be an 

exact l-o lower confidence point for e A . 

H is accepted iff F„ < F _ . An exact l-o lower 

o 2 — o:n 2 »n 3 

confidence point for $j/Ig 3 - 1/1 is 


3 1 „ °1 ~ °3 


So when F 2 < F^.^ ^ ue Kant 

s -1 

°1 a:n,,n, , 

L(F r F 2 ) ‘ JTT" 1-3 '' T 

This is not consistent with the firs, live coik.: ti<>.- 
the denominator of (2.5) to be I when condition (<} 
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Instead of using L(Fj,F 2 ) in equation (2.6), the simpler function 

In equation (2.7) will be used; which is obtained from equation (2.5) 

by setting b 0 = 0 . Note that Lim b, in equation (2.6) satisfies 
£ I— 2 

the requirement in condition 5). So the lower confidence point for e A is 


L A (F,,F 2 ) 


-i + r ‘ 


• f ) * ‘'-fal 


F 


_3_ 


-1 

1 


it s*/s*> 

^o:n^, 

( 

n 3 

• » if 

F a:n 1 , 

n 3 

Note that the confidence interval in equation (2.7) covers 

maximum likelihood estimator of 0^ 

, say 

e A . To demonstrate 

the M.L. estimator of , denoted by 9^ 

is 

e A - (F 1 -l)/(F 2 +n 1 ) if 

f, > i 


®A = 0 if 

F,il 


and the confidence interval is 



-1 + F F, + (1-F 

a : n^i a:n- 

F )F F 

\ *" a: V n 3 “ :n l ,n 3 1 . „ 

n, + F 

1 a:nj 

- F «», 

F "A 

* n 2 2 


1 
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- 1 * F -V' + l, - F «^ F « n , • * F > •' 


{1-F F „ ) < 0. 

0‘Dj»«° a.ni ^2 


"r 

r -i 


a: V 

S° L^F^.Fg) < 0 A < ® and @ A Is in the l-o confidence interval. 


1-a confidence interval for e A> given by 0 ± 0 A < u A ( r -j ,F 2 ) can be 
readily obtained by replacing o by l-o in equation (2.7). The 
upper l-o confidence point of the lower confidence interval is 
given by 


-1 + F 


-1 


U a (F,,F 2 ). 


l-a:n,,» 1 


F, + 0-F- 


-1 


> F Wn 


c-1 


U A (F r F 2 ) = 0 


"1 *3 

V n 3 


( 2 . 8 ) 


The procedure for determining this lower confidence point is 
similar to that described for the upper confidence point. 

From the formulas for upper and lower confidence intervals in 
equation (2.7) and (2.8), the two-sided confidence interval on r, 


1 
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can be readily obtained. The 1-o two-sided confidence interval for 
is given by LJ <_ o A <_ , where 

-1 + F I/9.n F 1 + 0" F /9.„ F / 9 .„ n ) F /9.„ „ 

a/Z:n^,« 1 a/Z:ni,« a/Zin^.n^ 7 a/Zrn^.n^ 1 


-1 + F l-o/Z:n r - F l + tl ” F l-a/Z:n 1 ,- F l-o/2:n 1 ,r. 3 )F l-a/2:n ] ,n 3 F l 
n 1 + F l-a/2:n,.- F 1-a/2:n,,n F 2 


This will be proved in Chapter 3. 

If either limit is negative it is replaced by iero. 

As mentioned in Chapter 1, due to the symmetry in the model, if 
a confidence interval fcr Is obtained, the confidence interval 
for e T can be easily obtained. The formulas for the confidence 
Intervals for e T are given below where Oj = Io t 2 /(o a 2 +o e 2 )0. The 
1-a lower confidence point for e T is 

•' * F .lL- F 2 * <'- F :L- F . : „„n > F ,:„,.n F i' 


The 1-a upper confidence point for e_ is 


-1 + F l-a:n 2 ,« F 2 + ^" F l-a:n 2 ,» F l-a:n 2 ,n 1 ^ F l-a:n ? ,n^ F 2 

VW- ~ r i -7^-—- <*•"» 

"9 n M 




The l-o two-sided confidence points for Oj are 

-1 + F o/2:n 2 ,- F 2 + (1 " F a/2:n k> - F a/2:n 2 ,n 3 )F a/2:n 2 ,n 3 F 2 1 

.L + F «/2:n 2 ,- F a/2:n 2 ,n 1 F l (2>12) 

+ F 1-l/2:n 2 ,- F 2 + (1-F l-a/ 2 :n ? ,- F l-«/ 2 :n 2 ,n 3 )F 1 -a/ 2 :n 2 ,n 3 F 2 

n 2 + F l-l/2:n 2 ,» F l-oi/2:n 2 ,n 1 F l 

If any.limit is negative it Is replaced with zero. 

From the confidence intervals obtained for o A and o T » and by 

equation (2.1) we get the corresponding confidence Intervals for p a 


The 1-a lower confidence interval for p a Is 


where = L a (Fj,F 2 ) wh ’ ch * s defined in equation (2.7). 
The 1-a upper confidence intervals for p a is 


( IU A \ 


where U A = u a^ F 1» F 2^ which is defined 1n e( l uation (2.8). 
The 1-a two-sided confidence interval for p A is 


(H \ / H 

[mwj - P A i [TTjq 


The 1-a lower confidence interval for p T is 
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where Lj = L T (,F 2 ) which Is defined in equation (2.10). 
The 1-a upper confidence interval for p T is 



•e Uj * which Is defined in equation (2.11). 

The 1-a two-sided confidence interval for 


(j * TUf) i p t i(rr-%) 


Where l)^ and are the confidence points given in the equation 

( 2 . 12 ). 
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3. Evaluation of the Methods 


In the preceding section confidence point(s) were determined 
for upper, lower and two-sided confidence intervals which have exact an 
confidence coefficient l-o for specified restrictions, and for general 
conditions the confidence coefficient should be close to l-o . In 
this section methods will be derived to determine how close 
the approximate confidence coefficients are to the specified ones. 

First consider the upper confidence Interval 

p C L (F r F 2 ) 1 ° A < -3 ‘ (3-D 

From equation (2.7) 


L(F 1 ,F 2 ) 


-1+F .1 Fn+d-F .I F )F „ F, 1 

o.-n,,- 1 o:n,,» ain^n,' o:n,,n- 1 

-J-L--L_3- 1 —£- if F > F 


1 o:n,,® o:n,,n 9 2 


n l *3 


If F, < F n 
1 o.n r n 3 

(3.2) 


From Section 1 

F i B s i 2 / s 3 2 for 1 * 1.2 ; also n i S 1 2 /e i * ^ 

for i = 1,2,3 are Independent chi-square random variables with n^ 
degrees of freedom respectively. 

Consider the following events 


A = {l(F r F 2 ) < e A ; F ] >0; F 2 > 0} 

B 1 - 1F 1 * F a: n ,.n 3 ' '2 
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B 2* {F l± F a: V n 3 l F 2 > °* L ( F l’ F 2>i°> (3.3) 

So B 2 B 0 and B| U « n where 

0 - ^ r F 2 , L ( F r F 2 )) : F 7 > 0; F ? > 0; L(F, ,F 2 ) > oj . 
Clearly P^(e^) * P(A). We will write events B 1 and B 2 In (3.3) - 


“1 - (F 1 > F a: n ,.n 3 > > B 2 ’ ' F 1 ± F .:n,< 3 -«> 


Prom this we obtain 


P L (e A ) = P(A|S 1 )P(B 1 ) * P(A|B 2 )P(B 2 ) 

W - p [L(f,.F 2 ) < 8 a |F, . F o:n)>n 3 ]P[F, > F a:n)i „ 3 ] 
♦P[L(F,.F 2 )<e A |F 1 <F o:n)ir , 3 ]P[F, iF, ;v „ 3 ) 


,tF «:p,.. F . : n 1 ,H3 )F g : n 1 ,n 3 F i 


1 o:n,a:n,,n 9 2 


*— e Ar l" am, ,n. 


x P[Fi > F _ ] + P[0 < e« IF, < F n ]P[F, < F n ] 
1 a.n^.n^- 3 *■ — A 1 1 — a:n^,n 3 J *■ 1 — a.nj.nj- 1 

But Is always greater or equal to zero so 



















1 


( Pi^ ^a:n i,n, IL^ 

(F “-r-' F -,.n 3 ) -^T- - 07 - F .:n,,.(’ + ",»A) 

PjUjVjy^ _ I 1 

n l n 3 / 0 A°1 P 2 F a: ni ,„ 2 iU 2| F ' >F ^",." 3 


p, . p _ Vl „ t (F «="i-~ F «="i."3 )f »:ni.n 3 

' " 1 V 2 F c,:n,,n 2 1 ’Vl p 2 e A F o:n,,n 2 U 1 

< ltn l°A> F a:r,,. | ] 

' "m F <,:„,,„, Us 1 Uj P ’ F « : ", ,n 3 


Let C, • (n 2 »,)/(n,e A » 2 


C 2 " t F «:fi|,~" F «:n 1 ,nj) F a:n,,n^( n 3 p l 9 2 9 A F o:n )1 n 2 ^ ^ 3 ' 6 * 
C 3 = ( ’ +n l e A )F o:n 1 ,/ (n li>2 e A F a:n ) ,n 2 ) 

Note C-j > 0 , C 2 < 0 and C 3 > 0. Then 
U 3 2 

P 1 • PC C 1 U , * C 2 -D7- - C3U3 < U 2 | f , > F a;n)i n 3 ] 

[WV U 3 > c 2 u 3 2 (u T *u 3 ) 2 c 3 u 3 

' ‘vv 1^7 °i ' W ‘ 1 3> 1 


r 







Let X = U-j + Uj . Y * , so (1-Y) = 


Then 


P, - P[C,YX + Cjd-YjV’x - C 3 (1-Y)X < U 2 |F 1 > F^^] 


p L (e A ) = Pf^YX + C 2 O-Y) 2 Y"'x - C 3 (1-Y)X<U 2 ; F, >f a:n ^ n l 
+ P[F, iF^^] (3.7) 


The random variables X, Y, Ug in (3.7) are independent and 

X~ x 2 (n 1 +n 2 ) , Y ~ b( j , U 2 ~ x 2 (n 2 ) 

Now consider the event {F, > F „ „ The following events are 
1 a:n,,n~ 3 


equivalent: 
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U + n 2»1 


n 2 p l 


Y > 


Fa:n T n 3 _ 

F a:n r n 3 + n 2 p l 


(3.8) 


Then (3.7) can be written as 

,P L (e A ) = [( c i Y + “ c 3 0-y)|x < u 2 ; 

+ p[y < F „ /(F „ +n,p,)l 

L — oin-jjny a:n-|,n 3 2 1 J 

= P f((C 1 +C 2 +C 3 )Y+C 2 Y~ 1 - (2C 2 +C 3 ))x < U 2 ; Y > a] + P[Y < A] 


= P[g(Y)X < U 2 ; Y > A] + P[Y < A] (3.9) 

where 

9<Y) = (C 1 +C 2 +C 3 )Y + C^' 1 - (2C 2 +C 3 ) 

With probability one Y e (0,1]. 

Consider the function 

(C^+Cjy 2 - (2C,+CJy + C, 
g(y)- - ----—-- 0 < y < 1 


g(y) = o iff (C 1 +C 2 +C 3 )y 2 - (2C 2 +C 3 )y + C 2 = 0; 


Clearly 




hence, g(y) has only two zeros for y e R, and they are 


. < 2 V C 3> * r C 3 2 - 4C ! C 2 ]V2 I* 

.1 " ~ 2(0^2+03) , j 

(2yc 3 ) - [C 3 ; -4C,C 2 ] 1 ^ I 

°2 " 2(C 1 +C 2 +C 3 ) f 

Also Lim g(y) - -® and Limg(y) = C, >0. Since g(y) is continuous 
y-*0 + y-*l 1 

in (0,1] then at least one of the roots (6 j, 6 2 ) 1n (0,1)* But 

= (C^C^) - C^" 2 = C 1 +C 3 +C 2 (l-y" 2 ) , 0 < y < 1 

But for y c (0,1), ^1 - < 0 . Also, (C^) >0 and ^<0, 

hence, if ye (0,1) then 

> 0 . 
dy 

So g(y) is monotonic increasing in (0,1]; also since 

Lim g(y) = - » and Lim g(y) = C, > 0 , only one of the two roots 
y-*0 + y>l 

can be in (0,1]. This gives the following two possible forms for g(y). 



Figure 3.1 







Figure 3.2 
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we obtain 

* P[g(Y)X < U 2 ; Y > A|A < Y < 6] P[A < Y < 6] 

+ P[g(Y)X < U 2 ; Y > A|Y > 6] P[Y > 6] 

' 1 • P[A < Y < 6] + [g(Y)X < U 2 ; Y > Ai Y > 6] 

“ P[A < Y < 6] + P[g(Y)X < U 2 ; Y > 6] 

- P y (6) - F Y< a ) + P ± u 2 ; Y ' 

where ' 

F Y (y) - P[Y < y] 

From the original probability in (3.9) we get 

P L (0 A ) = F y (6) - F y (A) + P[g(Y)X < U 2 i Y > 6] + F y (a) 

c F y ( 6) + P[g(Y)X < U 2 ; Y > «] (3.11) 

Also 

P[g(Y)X < U 2 ; Y > 6] = P[g(Y)X + X < U 2 + X; Y > 6] 

r U ? ♦ X 1 

• P[g(Y) + 1 < - ; Y > fij 

* p[[g(v) + l]" 1 ' Y > 5 ] (3. 12 ) 

Let Z = X/(X+Up) ; It follows that Z is distributed as Beta with 

/ n l +n 3 n 2 \ 

parameters - , -j-) • Also, let n^ be an even integer so that 

n^/2 = nij for 1 = 1,2,3 and hence the are positive integers. 

If X ] , X 2 , X^-re independent chi-square random 
variables with , m 2 , m^ degrees of freedom, respectively, and 
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then Yj, Y 2 are mutually independent Beta random variables with 
parameters (ny 2, ny2) and ((m 1 +m 2 )/2,m 3 /2^ respectively. 

From the result in (3.12) we obtain (let v denote [g(y) +1]" 1 ) 

p l (8 a ) * pftgm+i]' 1 > Z; Y > «] + f y («) 


/ 1 /* V B” 1 (m, ,m,) m, -1 itu-I m^nu-l m..-l 

J > 1 <'■*> 3 21 3 <’- 2 > 3 02 dy 

6 0 1 6 l 




(m,,m.) m,-l m 0 -l m,+m,-l 

■ v*> *j j nvM?* lUy) z 
6 0 


nyl 

l S 




m--r-l 

(l) r (- 2 ) 2 dz dy 


mi i m_-r-l 

v 1 B' 1 (m 1 ,m 3 )r(m 2 )(-1) 2 

* F Y (6) + r ? Q B(m 1 +m 3 ,m 2 )r(r+1) r (nyrf 

, btyj+i]' 1 

x // Vl (l-y)" 3 ' 1 / z <VW r - ? > dz dy 

6 0 

Vl B' 1 (m r m 3 )r(m 2 )(-l) m? r_1 
P L ( 0 A ) - F y(<5) + Ej B(ni 1 +m 3> m 2 )r(r+ fJrTnyr) 


i. 
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m,-l m,-l 

O-y) 3 


7 ’" 


r-XatyM] , 2 3 „y 

(n^+m^+tr^-r-l) ujr 


m 0 -1 


F («) ♦ T <-1) (rn^y 3 - 

F »' s) fit ispn 


,-1)1 


Tm^TJl/ Jr! (m^+ni^+m^-r-l) 


f' ;r’ (l-y) 3 

y (m.,+in ? +m,-r-l) dy 

« [g(y) + i] 1 2 3 


(3.13) 


From the expression In (3.13) w“ will compute P L ( 0 A ) 


Next consider the lower confidence Interval 
i»te A < U(F,.F 2 )] - P 0 (e A ) 
From equation (2.8) It follows that 


U(F r F 2 ) 


-i + f; 1 f, + (.i-f: 1 „ f, 

l-a:np» 1 v l-o:n^,» 1 -a:n^,i 


(3.14) 


»V l-a:n,,n, 1 


n, + F, F, F~ 

1 l-a:nj,«° 1-a.n^ ,n 2 <■ 


U(F r F 2 ) = 0 


” SlV * F l-«l.,.n, 


(3.15) 


Substitute equation (3.15) Into (3.14) to obtain 

Py(« A ) ■ p[[ 9(Y) ♦ l]' 1 I 2; Y > fij (3.16) 

where Y ~ Beta (m^rn-j) , and Z ~ Beta (mj+m^.n^); Y and Z 
are Independent. 









9(y) * (C 1 +C 2 +C 3 )y+ Cgy -1 - (C 3 +2C 2 ) 

6 Is the largest root of g(y). 

The C's in g(y) are obtained by changing the a's 
1-a's in the equation (3.6). So we have 

C Vi _ 

1 "lV2 F l-.:n 1 ,n 2 

^ F l-a:n 1 .-"^l-a:^ .n^l-arn^ ,n 3 

C2 n 3 p l p 2 9 A F l-o:n 1 ,n 2 

_< 1+n lV F l-«:n r . 

° 3 ”l' p 2 0 A F l-a:n 1 " i n 2 '" 

Note Ci > 0 , C 3 > 0; and C 2 > 0 if n^ > 2, n 2 > 

a < .10. 

Consider the function g(y) In equation (3.17) . Cl 
9(y) * o Iff (C,+C 2 +C 3 )y - ( 2 C 2 +C 3 )y + C 2 
hence, g(y) has only two zeros for y c R and they are 

. <SV C 3> * [c 3 2 -‘": 1 c 2 ] 1/2 
*i -rcy ♦ c 2 +t 3 j- 

(2C 2 .C 3 ) - tC 3 2 -4C,C 2 ]' /2 
‘i " 2(C, t C 2 '*Xp 


Also Urn g(y) 


and Lim g(y) ■ C, > 0 . 
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Since g(y) Is continuous in (0,1], then either both of the two roots 
of g(y) are in (0,1] , or neither of the roots are in the interval 
(0,1] , or the two roots are complex and not real. But A c (0,1] so 
that g(A) < 0 . For 


A ■ 



V 

n 3 


n _3 _ 

* °l n 2 


9(i) 


n*ir i r« 

A1 12 l-gin,,. l-g:n 1 ,n 3 
n 3 p l p 2 e A J? l-a:n 1 .n 2 


(3.19) 


0 


This follows from the proof of (b) In Equation (3.10). 

From the above it Is clear that the two roots of g(y) are in 
(0,1], which results in the following form for g(y). 



From equation (3.16) we evaluate Py(e A ) as follows 
_ , m,-r-l 

V' (-1) 2 (i^+iyuyl)! 

P U^ e A^ * (m^-l)! (rn^-1)! (m^-r-l) !rT(m^^^+m^r-TT 


L 
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The two-sided confidence interval Is 

PCL 1 (F 1 ,F 2 ) < 0 A < U 1 (F 1 ,F 2 )3 = P 2 (e A ) (3.21) 

where (F^*Fg) is given by (3.2) with a replaced by a/2 and 
U|(F].F^) given by (3.15) with l-o replaced by l-o/2 . 

To show (3.21) let A and B be the following events 

A ■ {L-j(F] .Fg) £ 6^) B “ (0^ £ Uj(Fj *^ 2 ^ 

Then 

P 2 (0 A ) * P(Af\B) (3-22) 

From (3.22) we obtain 

P 2 (e A> = P(A( ^ B ) * 1 - P(ATiB) - 1 - P(Al)B) 

= 1 - [P(*) + P(B) - P(AflB)) 

• 1 - P(*) - P(B) ♦ P(Knff) (3.23) 

Also 

P<*) * 1 - P L (e A ). P(B) - 1 - P u (e A ) 
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Lemma 3 . F ]_ o;n >n is Increasing in n 2 if n^ > 2 
n, > 2, a < .10 . 


F a:nf,n 2 < F a:n r n 2 1f n T > n i 0>y Lemma 2). 


F' 1 < F' 1 

l-a:n 2 ,n*- l-o 

And from this the result follows. 

Using the Lemmas in expression (3.24) we have 
1) C 2 < Cj since 


F l-a:n 1 ,. F l-a:n 1 ,n 2 i F a:n 1 .- F a:n 1 ,n 


- 1 + AiF, + Bt F7 1 - 1 + A,F, + B, 

c , i ’ • c 2 i“ -npnfc 11 1 


2) f l:n v i F i-a:n,, s ° A 1 i A 2 


-1 ♦ A,F, + BjFj 1 - 1 * A Z F 1 ♦ BjFj' 
"l * ~ "2 + C 2 F 2 

3) F a:n,,- F «:n r n 3 - F l-a.-n, ,- F l-a:n, .n. 


a:n-|,» ain^.n^ - l-a:n^,» l-a:n^,n 3 


(1 " F »:n 1 ,- F a:n 1 .n 3 ) - (1_F i-a:n 1 ^l-a:^ ,n 3 
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,«^a:nj ,n,^o:n^ .n^ — ^ l-ain^.n^ 

then B, < 


-' + ¥-, * ¥1* , -1 ♦ Vi ♦ 6 2 F i 1 

’ "l + C 2 f 2 


"l + C 1 F 2 
Going back to (3.23) we get 

P. 


2 <V * 1 * ( 1 - p L (e A ) ) - KM * 0 
P 2 (e A ) = P L (e A ) + ?„(«*) - l 


(3.25) 


Numerical Integration will be used to evaluate the Integrals 
In P L ( 0 A ) , P u( 0 /\)* and P 2^ e A^* The method used is Simpson's Rule 
which gives 


J ^ f(y)<iy~| 


f(a) + 4f(a+h) + 2f(a+2h) 


+ 4f(a+3h) + 2f(a+4n) + ... + 2f(a+(n-2)h) + 4f(a+(n-l)h) 

♦ f(D)] 

where h Is the length of the subintervals in which (a,b) has been 
partitioned. We used h = .001. 

Various combinations of sample sizes were cnosen, and for each 
of the combinations several sets of , and p 2 , and different 
o's were used. 
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A list of the values of I, J, p^, p^, and a used Is given 

below 

(I.J) : (3,3);(3,7);(3,11);(5,3);(5,5);(5,11);(7,3); 

(7»7);(7,11);(9,3);(9,9);(9,11);(11,3);(11,7);(11,11) 
o A ? /o £ 2 : .003, 1, 10, 100 (4.1) 

o T 2 /a e 2 : 0, 1, 10, 100 
1-a : .90, .95, .975, .99 
d } ■ 1 + o a A 2 /o e 2 

P 2 » 1 + I o T 2 /o e 2 

The results are given in Tables 2, 3, and 4. The first five columns 
in each table show the I, J, n^, n 2 , and n^ values. Each of the remaining 
columns give the results for the specified 1-a confidence coefficients. 
For each of these columns the ranges that the confidence coefficients vary 
(as p 1 and p 2 vary over their allowable values) are shown. 
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The approximate confidence coefficients corresponding to the 
lower confidence intervals are very close to the specified ones, shown 
in Table 2. The approximate confidence coefficients for the upper 
confidence intervals range from values very close to the specified ones 
to some high values, shown in Table 3. Similar results hold for the 
two-sided confidence limits shown in Table 4. 

For I fixed the range of the confidence coefficients is closer to 
the nominal value when J increases. For J fixed the upper bound of 
the range is higher when I increases. So the confidence coefficient 
will be better when J is large. However, in general the approximate 
confidence procedures seem adequate if I and J are each 5 or bigger. 




Ranges of confidence coefficients for lower confidence intervals on 
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Table 4 

Range of confidence coefficients for a two-sided confidence 
interval on 


I 

J 

n l 

n 2 

n 3 

1 - a 

= .90 

1 - a 

= .95 

3 

3 

2 

O 

4 

.90103 - 

.93208 

.95072 - 

. .97000 

3 

7 

2 

6 

12 

.90043 - 

.91349 

.95030 - 

.96000 

3 

11 

2 

10 

20 

.89950 - 

.90849 

.94827 - 

.95702 

5 

3 

4 

2 

8 

.90096 - 

.93896 

.95062 - 

.97273 

5 , 

5 

4 

4 

16 

.90077 - 

.92684 

.95051 - 

.96725 

5 

11 

4 

10 

40 

.89929 - 

.91272 

.95031 - 

.95884 

7 

3 

6 

2 

12 

.90094 - 

.94310 

.95604 - 

.97436 

7 

7 

6 

6 

36 

.89871 - 

.92475 

.94984 - 

.96607 

7 

11 

6 

10 

60 

.89896 - 

.91673 

.94955 - 

.96130 

9 

3 

8 

2 

16 

.90092 - 

.94619 

.95059 - 

.97647 

9 

9 

8 

8 

64 

.89835 - 

.92374 

.94929 - 

.96546 

9 

11 

8 

10 

80 

.89878 - 

.92016 

.94938 - 

.96334 

11 

3 

10 

2 

20 

.90090 - 

.94863 

.95058 - 

.97782 

11 

7 

10 

6 

60 

.89790 - 

.93203 

.94926 - 

.97006 

11 

11 

10 

10 

100 

.89868 - 

.92315 

.94927 - 

.96510 
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